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Arguably the most intriguing aspect of the 
physics of cuprates is the close proximity between 
the record high- Tc superconductivity (HTSC) and 
the antiferromagnetic charge-transfer insulating 
state driven by Mott-like electron correlations. 
These are responsible for the intimate connec- 
tion between high and low-energy scale physics 
[l-3j, and their key role in the mechanism of 
HTSC was conjectured very early on [4j. More 
recently, the detection of quantum oscillations in 
high-magnetic field experiments on YBa2Cu306+x 
(YBCO) has suggested the existence of a Fermi 
surface of well-defined quasiparticles in under- 
doped cuprates [5, 6j, lending support to the 
alternative proposal that HTSC might emerge 
from a Fermi liquid across the whole cuprate 
phase diagram [7, 8j. Discriminating between 
these orthogonal scenarios hinges on the quan- 
titative determination of the elusive quasiparti- 
cle weight Z, over a wide range of hole-doping 
p. By means of angle-resolved photoemission 
spectroscopy (ARPES) on in situ doped YBCO 
[9j, and following the evolution of bilayer band- 
splitting, we show that the overdoped metal elec- 
tronic structure (0.25 <p< 0.37) is in remarkable 
agreement with density functional theory [10-12J 
and the Z = 2p/{p-\-l) mean-field prediction [13, 14j. 
Below p^ 0.10-0. 15, we observe the vanishing of the 
nodal quasiparticle weight Z^r; this marks a clear 
departure from Fermi liquid behaviour and — con- 
sistent with dynamical mean-field theory [15| — is 
even a more rapid crossover to the Mott physics 
than expected for the doped resonating valence 
bond (RVB) spin liquid fT3llT4] . 

Formally, the degree of quasiparticle integrity is re- 
vealed by Z/c = jAcoh{k^oj) duj^ i.e. the integrated spec- 
tral weight of the coherent part of the single-particle 
spectral function A{k,uj) = Acoh{k^oj) -\- Aincoh{k^oj) 
probed by ARPES [Ml [17]. Experimentally, while in the 
optimally-to-overdoped regime Zk is believed to be finite 
- yet quantitatively undetermined - at all momenta both 
above and below Tc^ the situation is much more contro- 
versial in the optimally-to-underdoped regime [I6l [17]. 
Although it has been conjectured that the T = ex- 
trapolation of the pseudogap state is a nodal (N) liquid 
[m [E], a determination of Zn has not been possible. 
As for the antinodal (AN) quasiparticle spectral weight. 



the ARPES spectra are characterized by a dramatic tem- 
perature dependence, with broad incoherent features in 
the pseudogap state and quasiparticle-like excitations 
emerging below Tc [iSl [17] ; this observation led to the 
early proposal that the onset of HTSC might be thought 
of as a 'coherence transition' [20 . Detailed doping- 
dependent studies reported a decreasing Zan upon re- 
ducing p [2T | [22 ] : however, whether Zan^P [22 or van- 
ishes at finite p [21] remains unresolved, because of the 
experimental difficulty in the quantitative discrimination 
between Acoh and Aincoh- As discussed in greater de- 
tail in the Supplementary Information, where we present 
an analysis of the quasiparticle spectral weight across the 
whole YBCO phase diagram, this is particularly challeng- 
ing on the underdoped side where Zj^ vanishes; and even 
in the overdoped regime, this allows at best a relative - 
rather than absolute - determination of Z^. 

As an alternative, potentially more quantitative ap- 
proach, in materials with Cu02 bilayers within the unit 
ceh such as Bi2Sr2CaCu208+5 (Bi2212) and YBCO, 
the quasiparticle strength Zk might be estimated from 
the bonding (B) and antibonding (AB) band splitting 
e^'^^{k) = e{k) ^ t±{k), where e{k) is the quasiparti- 
cle dispersion with respect to the in-plane momentum /c, 
and t±{k) accounts for the interplane coupling. Access 
to Zk is provided by the fact that both e{k) and t_\_{k)^ as 
measured by ARPES, correspond to effective quantities; 
under the assumption that correlation effects renormal- 
ize the bandwidth but not the functional form of the 
quasiparticle dispersion, these effective quantities are re- 
lated to the density functional theory values through the 
same renormalized quasiparticle strength Zk. This as- 
sumption is most appropriate along the nodal direction 
of cuprates [I6l [17] , which is the focus of our study since 
nodal quasiparticles - being unaffected by the opening of 
pseudo and superconducting gaps - provide access to the 
normal state underlying the emergence of HTSC. This 
approach for the extraction of Zk will be experimentally 
validated based on a direct comparison with results for 
the quasiparticle spectral weight. 

The renormalization of the in-plane hopping probabil- 
ity due to correlations is given by the ratio between the 
number of ways an electron can be added to the Mott-like 
system as compared to the independent particle system. 
While in the latter this is proportional to the total num- 
ber l-\-p of available band states [1 at half-filling (p = 0), 
and 2 for a completely empty band (p = l)], in the Mott- 
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FIG. 1: Fermi surface and band dispersion across the 
YBCO phase diagram, (a) Schematics of the Fermi surface 
of optimahy doped YBCO from band-structure calculations: 
three sheets of Fermi surface {FSch, FSb, FSab) are de- 
rived from the one-dimensional CuO-chain (Ch) band, and the 
bonding (B) and antibonding (AB) Cu02-plane bands, (b) 
YBCO doping-range for bulk (oxygen content, top axis) and 
surface (hole-doping p, bottom axis) investigated in this work. 
Color dots refer to different samples of various oxygen con- 
tent; from the combination of as-cleaved and K-deposited sur- 
faces, a doping p ranging from 0.37 to 0.02 could be accessed, 
(c-h) Fermi surface and band dispersion, as determined by 
ARPES at T = 20K, for three representative hole-doping lev- 
els obtained on: self-doped as-cleaved YBC07 [(c,d), p = 0.29, 
red dot in b], a second K-evaporated YBC07 sample [(e,f), 
p = 0.20, purple dot in b], and heavily K-evaporated YBC06.3 
[(g,h), p = 0.06, blue dot in b]. At T = 20K the data are 
representative of the overdoped normal metal (c, large FSb 
FSab barrels), the nearly optimally doped superconductor 
(e, gapped Fermi surface with nodal Dirac points), and the 
underdoped pseudogap state (g, Fermi arcs). 



correlated case it is proportional to the doping p away 
from half-filling, times a factor of 2 because for each hole 
induced by doping an electron can be added with spin 
up or spin down [1 . This leads to the mean- field result 



Z = 2p/{p^l) [131 [14], aiid is the same for in-plane and 
interplane hopping since both processes are governed by 
the same correlation effects. The effective bilayer split- 
ting can be written as Aef '^^ = 2t^{k) = 2Zkt^^^{k); 
its detection by ARPES and the a-priori knowledge of 
t±^^{k)j as derived from local-density approximation 
(LDA) band structure calculations, allows a quantitative 
determination of Zfc. 

A study of this kind has been performed on Bi2212 at 
the antinodes in the optimal-to-overdoped regime [23] : 
however, a systematic investigation in the underdoped 
regime is lacking. A potentially more promising candi- 
date for such a study of Z^ is YBCO, which is charac- 
terized by a larger bilayer splitting especially along the 
nodal direction [TQtfl2] : furthermore, although YBCO 
has long been considered unsuitable for ARPES because 
of the polar-catastrophe-driven overdoping of the cleaved 
surface [Q] |12] , the in situ doping method based on the 
deposition of potassium allows exploration of a doping 
range much wider than for any other cuprate family [9|. 
Here we apply this approach to detwinned YBa2Cu306+aj 
single crystals for several values of the bulk oxygen con- 
tent: (6+x) = 6.34, 6.35, 6.51, and 6.99 (hereafter referred 
to as YBC06.3, YBC06.5, and YBC07). As summa- 
rized in the YBCO phase diagram of Fig.[T]3, by perform- 
ing ARPES on different as-cleaved and K-deposited sam- 
ples we can follow the evolution of the electronic structure 
from the heavily overdoped (p:^0.37, by far the highest 
among all overdoped cuprates including Tl2Ba2Cu06+(5 
[24]), to the deeply underdoped regime {pc^0.02). 

In the overdoped regime (p = 0.29, Fig.jl]^), we observe 
three Fermi surfaces {FSch^ FSb, FSab), in agreement 
with LDA band-structure calculations (Fig. [Ik). These 
correspond to the one-dimensional CuO-chain (Ch) band, 
and the B and AB Cu02-plane bands. The dispersions 
of the bands seen by ARPES (Fig.jlji) are also consistent 
with the LDA results, with the exception of the BaO 
band that in the calculations is located at the Fermi en- 
ergy [TOldT], but in the experiment is found at ~450 meV 
binding energy (similar disagreement is encountered for 
TIO t24| and BiO bands [25^ in Tl- and Bi-cuprates). 
Upon reaching the underdoped regime by potassium de- 
position on the as-cleaved surfaces, behaviour consis- 
tent with the established phenomenology of underdoped 
cuprates is seen [TS] [26H29] : the Cu02-plane Fermi sur- 
face barrels become partially gapped and reduce to four 
nodal 'Dirac points' (somewhat broadened by resolution) 
in the p = 0.20 superconductor (Fig.[l^), and to four 
extended Fermi arcs in the p = 0.06 pseudogap phase 
(Fig. [It). The spectral function becomes progressively 
more incoherent, with reduced low-energy quasiparticle 
intensity and a continuum of spectral weight extending 
to increasingly higher binding energy (Fig.[T]F,h). We also 
observe a progressive change in FSch volume and BaO 
binding energy (Fig.[Tfc-h); as discussed in the Supple- 
mentary Information, monitoring the evolution of these 
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FIG. 2: Doping evolution of the gap in YBCO. (a-c) 

ARPES image plots and selected energy distribution curves 
(EDCs), from the antinodal to the nodal region, along the 
AB-band Fermi surface and minimum gap contour. The data 
were obtained at T = 20K on the same three samples of 
Fig.[l] (a) overdoped normal metal {p = 0.29, T > Tc); (b) 
nearly-optimally doped superconductor (p = 0.20, T < Tc); 
(c) underdoped pseudogap state {p = 0.06, Tc < T < T*). 
While gapless metallic behavior is observed for p = 0.29 (a), 
a v-shape gap consistent with the d-wave functional form is 
detected for p = 0.20 (d), as evidenced by the evolution of 
both superconducting peak (SCP) and leading-edge midpoint 
(LEM) as a function of the Fermi surface angle 9 defined in 
Fig.TJi. For p = 0.06 (e), a much more rounded gap profile 
is observed, consistent with the established pseudogap and 
Fermi arc phenomenology (in this latter case SCP corresponds 
to the minima in the second energy derivative of the EDCs). 



features allows the precise determination of the effective 
surface hole-doping p induced by potassium deposition. 

As for the superconductivity related features, Fig.|2] 
shows image plots and selected energy distribution curves 
(EDCs) from the antinodal to the nodal region, along the 
AB-band Fermi surface and minimum gap contour. The 
p = 0.29 EDCs have no marked momentum dependence, 
typical of metallic behavior (Fig.|2|L). At p = 0.20 (Fig.|2]o 
and d), a clear v-shape gap consistent with the dx2-y2 
HTSC gap is evidenced by the momentum evolution of 
both leading-edge midpoint (LEM) and superconducting 
peak (SCP). The A | cos {k x) — cos (ky) \ fit of the average 



of the two (Fig.|2]i) returns a maximum antinodal gap 
value A c^ 34 meV, consistent with the 2 A ~ 70 meV re- 
ported for HTSCs with T^^^^ - 95 K [29]; this also gives 
^f/^2 = 15 ±1, where V2 is the slope of the gap at the 
node and vf = 1.42 ± 0.05 eVA. At p = 0.06, the quasi- 
particle peaks and d-wave-like superconducting gap are 
replaced by broad incoherent features with marginally 
defined LEM and SCP (Fig.|2]3); these describe a rounded 
pseudo-gap-like profile [27], with a much larger A at the 
antinodes (Fig.[2fe). This behavior, and the associated 
Fermi arc phenomenology (Fig. fit), are hallmarks of the 
pseudogap phase [18l [26H291 . 

The (T = 20 K) results in Fig.[l] and |2] reproduce the 
major features of the cuprate electronic structure in the 
overdoped normal metal {p = 0.29, T > Tc), the nearly- 
optimally doped superconductor (p = 0.20, T<Tc), and 
the underdoped pseudogap state (p = 0.06, Tc<T <T^). 
Most importantly, they demonstrate that in situ potas- 
sium doping of cleaved YBCO provides a very effec- 
tive means of exploring the whole cuprate phase dia- 
gram via surface-sensitive single-particle spectroscopies. 
We will now discuss the evolution of kp and B-AB bi- 
layer band splitting along the nodal direction, aiming 
at a quantitative determination of Zk. As shown in 
Fig.[3^ for YBC06.3, the B and AB bands are clearly 
resolved in the overdoped regime (p = 0.24), while only 
one dispersive feature is detected at low doping values 
(p = 0.06); a similar evolution has been observed on all 
YBCO samples. By fitting the momentum distribution 
curves (MDCs) at Ep with Lorentzian lineshapes con- 
voluted with a Gaussian to account for the instrumental 
resolution (Fig.l3k,b), we extract kp^N from p = 0.37 all 
the way to 0.02. Bilayer splitting is progressively re- 
duced upon underdoping and below p :^ 0.15 only one 
Cu02 band is detected (Fig.pp). From the linear fit of B 
and AB nodal kp^N versus doping we identify a vanishing 
Pc = 0.12 ±0.02, consistent with the recent interpretation 
of optical - and thus bulk sensitive - data from YBCO, 
which also suggest the suppression of bilayer splitting be- 
tween 6.5 and 6.75 oxygen content [30 . 

As for the evolution of Z/e, we start from the heav- 
ily overdoped regime {p > 0.25, region III in Fig.pp), 
where ARPES is in overall quantitative agreement with 
the LDA calculations. At p = 0.37 (i.e., overdoped 
normal metal), we can accurately determine both the 
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150 ± lOmeV. From the LDA values 
120 meV and t^^^{AN) c^ 150 meV (Sup- 
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plementary Information), and making use of the rela- 
tion Aef '^^ = 2Zkt^^^{k), we obtain Zn ^ 0.54 and 
Zan — 0.50. Upon underdoping, we focus on the nodal 
direction where the bilayer splitting can be detected both 
as A/c^'^ and Ae^' , owing to the presence of the gap 
nodes. While the results in region II are still at least qual- 
itatively consistent with LDA (0.12 <p < 0.25), the dis- 
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FIG. 3: Doping evolution of nodal kpj^, ^F,Ar? and Zn- (a) ARPES nodal dispersion and corresponding Ef momentum 
distribution curves (MDCs; ±15meV integration, see shaded region), from as-cleaved and K-deposited YBC06.34, for light 
polarization parallel to F-S. (b) Similar MDCs from YBC07, for light polarization parallel to F-Y (note the polarization- 
dependent change of Ch, B, and AB intensity), (c) Evolution of kpj^ (down triangles) and kpj^ (up triangles) along the 
nodal direction (distance from the F point in percentage of F-S), as determined by ARPES across the entire hole-doping 
phase diagram of YBCO; below p:^0.15 the bilayer B-AB splitting vanishes and only one single kF,N is detected (diamonds). 
(d) Nodal quasiparticle renormalization Zn as determined from the observed bilayer splitting (Zn = ^^n /'^t^^ (N) with 
t;5l^'^(A^)^120meV) and the low-energy spectral- weight ratio (SWR) defined as f~'^ I{kF,N,(jo) duj/ f~^^y I{kF,N,(jo) duj from 
the ARPES intensity I{k,(jj) = Io{k)f{(jj)A{k,(jj). The SWR was rescaled so that the two curves match in the 0.23—0.37 doping 
range (see Supplementary Information). Thick gray and orange lines are spline guides to the eye, while the dashed red line 
is the mean- field relation 2p/(p-\-l) [131 I14j. For the bilayer splitting data, the error bars are defined from the B-AB MDC 
fits when bilayer splitting is detected, and from the experimental resolutions otherwise; for the SWR data, they are calculated 
from the spread in the SWR values for integration windows smaller than [0.1 eV, — cxo] (Fig. 3 in Supplementary Information). 



appearance of bilayer splitting in the underdoped regime 
(p< 0.12, region I in Fig.lsfc) marks a clear departure from 
the independent-particle description. Note that LDA 



would give a nearly doping-independent A/c^'^ with 
(^F,Ar; ^f,Ar) = (33.0; 38.7), (35.2; 39.5), and (38.0;42.2) 
for 07. 0, 06.5, and 06. 0, respectively [11 . Furthermore, 
the value of bilayer splitting is also rather insensitive to 
electrostatic effects associated with the presence of differ- 
ent concentrations of 0^~ ions in the CuO chain layers or 
K^+ above the BaO terminations, due to the vanishing 
potential these charged layers exert on the Cu02 bilayer 
(Supplementary Information). 

The bilayer-splitting-derived results for Zn across the 
whole phase diagram are presented in Fig.[3]i, together 
with the nodal spectral- weight ratio (SWR) estimated 
at Ef (see caption for definition). When calculated for 
integration windows smaller than [O.leV, — oo] in bind- 
ing energy, the SWR provides a relative measure of the 
degree of quasiparticle integrity Zn (Supplementary In- 
formation). Although not a quantitative estimate of Zn, 
the doping evolution of the low-energy SWR is reminis- 
cent of the quasiparticle renormalization inferred from 
the bilayer splitting. When the two results are plot- 



ted together (after the necessary rescaling of the SWR, 
so that the two curves match in the 0.23 — 0.37 doping 
range), the similarity is striking (Fig.pd), validating the 
quantitative determination of Zn from the evolution of 
bilayer band splitting. Furthermore, the agreement with 
the 2p/{p-\-l) relation above p:^ 0.23 provides an unprece- 
dented quantitative estimate of quasiparticle integrity in 
overdoped cuprates, which supports a description of the 
normal-state electronic structure based on density func- 
tional theory with nearly isotropic quasiparticle renor- 
malization factors, as in the spirit of Fermi liquid theory. 
Below p c^ 0.23, however, Zn deviates from 2p/(p+l), 
vanishing in the p ^ 0.10 — 0.15 range where the error 
bars in Fig.[3]i define an upper limit Z ~ 0.065. From 
this we can estimate a corresponding upper limit for the 
nodal bilayer splitting Ae^' c^l5.6meV, which is con- 
sistent with the value of SmeV reported based on the 
analysis of quantum oscillation results [ST, "32^ . Note also 
that in region I, not only have the bilayer splitting and 
Zn vanished, the doping evolution of kp^N appears to 
have saturated, despite the fact that the Cu02 planes 
are still becoming more underdoped as evidenced by the 
further increase of pseudogap magnitude and high-energy 



incoherent spectral weight. This indicates that, upon ap- 
proaching the charge- transfer insulator, the quasiparticle 
integrity is lost even faster than hinted by renormalized 
mean-field treatments of the doped spin liquid |13l[T4]. It 
also suggests that the emergence of the Mott physics in 
cuprates might be captured by cluster dynamical mean- 
field theory calculations, which predict a similarly rapid 
suppression of Zn in the pc^ 0.10— 0.15 doping range [15 . 
Finally, we note that a vanishing Z^ does not im- 
ply vanishing of the total single-particle spectral weight, 
but rather a redistribution from coherent to incoher- 
ent components of A{k^uj). While for p = 0.06 we find 
Zn — 0, spectral weight is still present all the way to 
Ef (Fig.[3|i and Supplementary Information); however, 
the distinction between Acoh and Aincoh is becoming 
completely blurred. This suggests something similar to 
early exact diagonalization cluster calculations for Mott- 
Hubbard and charge-transfer models J], where it is the 
total and not just the coherent spectral weight that van- 
ishes as 2p/(p+l). Indeed, optical spectroscopy [33] 
shows that the total low- frequency {uj < 75meV) inte- 
grated spectral weight - whether associated with a well 
defined Drude peak or an incoherent band - scales to zero 
with the doping p [33] . More fundamentally, the loss of 
quasiparticle integrity implies the breakdown of concepts 
such as Fermi surface, band dispersion and Fermi veloc- 
ity, as well as Luttinger's sum rule, and suggests that be- 
low p<0.10 — 0.15 the physical properties of underdoped 
cuprates (in zero-field) are dominated by incoherent ex- 
citations. 



slum. We acquired data on the as-cleaved surfaces, as 
well as after subsequent potassium evaporations from 
a commercial SAES alkali- metal dispenser; evaporation 
current and time were precisely monitored (in the 6.5 
to 6.7 A and 30 to 120s range, respectively). The tem- 
perature was kept constant at 20 K for the duration of 
the experiments to maintain the most stable experimen- 
tal conditions before, during, and between subsequent 
potassium evaporations. 
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METHODS 

Sample preparation. For this study we used de- 
twinned YBa2Cu306+x single crystals with x = 0.34, 
0.35, 0.51, and 0.99 (see Ref. 9 for details on growth, an- 
nealing, and detwinning procedures). 
ARPES experiments. ARPES measurements were 
performed at Beamline 7.01 of the Advanced Light 
Source, with linearly-polarized 100 and 110 eV photons 
and a Scienta- R4000 electron analyzer {AE c^ 30meV 
and A^ c:^ 0.1°). The light polarization was parallel to 
the plane of emission, while the samples were aligned 
with the in-plane Cu-O bonds either parallel or at 45° 
with respect to the polarization (resulting in more in- 
tensity for the antinodal/chain or nodal features, respec- 
tively). Fermi surface maps were obtained by integrat- 
ing the ARPES intensity over a 30meV energy window 
about Ef^ and then normalized relative to one another 
for display purposes. Mounted on a 5-axis manipulator, 
the samples were cleaved at 20 K and 2.5x10"^^ torr, and 
oriented in situ by taking fast Fermi surface scans. 

The self-doping of the YBCO surfaces for various 
bulk oxygen contents was tuned towards the underdoped 
regime by in situ deposition of submonolayers of potas- 



SUPPLEMENTARY INFORMATION 

Effective hole-doping determination 

As discussed in the paper, the data presented in Fig. 1 
and 2 (paper), consistent with the observed doping and 
the T = 20K temperature of the experiments, repro- 
duce the major features of the cuprate electronic struc- 
ture in the overdoped normal metal {p = 0.29, T > Tc), 
the nearly-optimally doped superconductor {p = 0.20, 
T <Tc)^ and the underdoped pseudogap state (p = 0.06, 
Tc<T <T^). Most importantly, they demonstrate that 
in situ potassium doping of the cleaved surface of YBCO 
is a very effective mean of exploring the whole cuprate 
phase diagram via surface-sensitive single-particle spec- 
troscopies. To this end, a crucial aspect of the data 
analysis is the conversion from K-evaporation to the 
effective hole doping p per planar copper. If each K 
atom deposited on as-cleaved YBCO donates one elec- 
tron to the surface electronic structure, the resulting 
K^^ ions should be more likely located on the insulating 
and charge-neutral BaO termination, rather than on the 
(l+)-charged CuO chain termination. This is confirmed 



by our scanning tunneling microscopy (STM) study of 
various YBCO samples after K-evaporation, which con- 
sistently finds K atoms on BaO terminations only; this 
is shown in Fig.|4^ for the case of ortho-II YBC06.5. 
STM was performed at 10 K and 2.0x10"^^ torr with a 
beetle type STM, and the same SAES getter source and 
procedures used in the ARPES experiments; tungsten 
tips were cleaned in situ by Ar sputtering and sharpened 
through gentle contact with an Au surface. Topographic 
images on BaO (CuO chain) terminations were acquired 
with bias voltage and tunnelling current of —200 mV and 
200pA (-600 mV and 500 pA). We studied a total of 10 
samples of different bulk doping from the same batches 
as those used for ARPES. 

With respect to the doping evolution of the top-most 
BaO, CuO-chain, and Cu02-plane layers, the ansatz to 
be proven is that the electrons donated by potassium 
will be sufficiently delocalized to lead to a uniform car- 
rier doping. In Fig.[4]3 we show, for various as-cleaved 
and lightly K-deposited samples, the binding energy of 
the BaO band (BEBao)^ as well as the volume (counting 
electrons) of CuO-chain and (bonding and antibonding) 
Cu02-plane Fermi surfaces {FSch^ FSb^ and FSab)- 
Note that, in this very overdoped regime in which all 
Fermi surfaces are well-defined, the effective hole doping 
p can be determined directly from the average volume 
of FSb and FSabi with no ambiguity. Interestingly, 
the same quantities estimated from LDA calculations for 
fully oxygenated YBCO 7, upon a rigid shift of the chem- 
ical potential, show good overall quantitative agreement 
suggesting that, beyond :^0.25, the electronic structure is 
well described by independent-particle density-functional 
theory (for the BaO band, which LDA puts closer to 
the Fermi energy, the values in Fig.[4]3 were obtained by 
shifting it down self- consistently to match the BE ob- 
served for p = 0.37 as the extreme starting point). As 
we further hole-underdoped the system below p c^ 0.25 
by K-evaporation, we cannot rely on the volume of FSb 
and FSab in determining p because of the opening of 
the superconducting gap and, eventually, the collapse of 
the Cu02 Fermi surfaces into disconnected Fermi arcs. 
We can however follow very clearly the increase of the 
(electron) FSch volume (Fig.|4^ for YBC07, sample 1) 
and of the BaO-band binding energy [Fig. ld,f,h (paper), 
for various bulk dopings]. The results for all evapora- 
tion steps on YBC07S1 are summarized in Fig.|4]3, and 
demonstrate that both FSch and BEBaO scale linearly 
with the amount of deposited potassium, all the way to 
the highest surface doping. Similar scalings are observed 
for all bulk oxygen contents, leading to one universal lin- 
ear scaling relation between AFSch and ABEBaO^ for 
all YBCO samples (Fig.[4]i). Note that the increase in 
both FSch and BEBaO upon K-evaporation (i.e., up- 
ward slopes in Fig.[4fc,d) confirms once again the electron- 
nature of the doping induced by K (i.e., effective hole- 
underdoping) . Furthermore, the linearity over the whole 



doping range is consistent with the uncorrelated be- 
haviour predicted for both BaO and CuO chain-derived 
electronic states in LDA+U calculations [11]. Quantita- 
tively, the K-induced rate ABE^ao/AFS'c/i = 0.0323 ± 
0.0012 [eV(%5Z)-i] obtained from Fig.jiji matches al- 
most perfectly the ratio ~ 0.0331 obtained from the 
evolution of ABEBaO and AFSch versus hole doping 
p in both as-cleaved results and LDA calculations in 
Fig.[4]3 [here we observed AB EBaO /p = 0.b0b9 {eV p-'^) 
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FIG. 4: K-evaporation: effective hole-doping deter- 
mination, (a) Topographic STM images (4x4 nm^) from 
oxygen-ordered ortho-II YBC06.5 confirm the approximately 
50-50% CuO-BaO termination of the cleaved surface, the 
ortho-II alternation of full and empty chains, and the pref- 
erential location on the BaO termination of evaporated K 
(red circle), (b) Binding energy of the BaO band (BEBao), 
and volume (counting electrons) of FSch, FSb, and FSab, 
plotted versus the hole doping p together with the results of 
LDA band-structure calculations [TT] . In this very overdoped 
regime, in which all Fermi surfaces are well-defined on both as- 
cleaved and moderately K-deposited YBCO, the hole doping 
p is determined from the average volume of FSb and FSab 
(symbols refer to different bulk oxygen contents) . (c) Volume 
of FSch and variation of the BaO binding energy (ABEbuo) 
versus the K-evaporation time for YBC07. (d) The relative 
evolution of ABEBaO and AFSch upon K-deposition is lin- 
ear over the whole range for all samples, (e) Fermi surface 
mapping from as-cleaved and K-deposited YBC07. 



Structure/doping 


Nodal bilayer splitting (eV) 


Bulk YBCOe.o 


0.19 


Bulk YBC06.5 


0.21 


Bulk YBCO7.0 


0.24 


Slab YBCO7.0 




B&Osurf 


0.18 


CuOsurf 


0.23 


Slab YBCO7.0 with K 




BSiOsurf 


0.21 


CuOsurf 


0.19 



TABLE I: Nodal bilayer splitting values for various YBCO 
bulk and surface structures, as obtained from WIEN2K den- 
sity functional theory calculations for kz=0. 



and AFSch/p = 15.269 {%BZ p'^)]. This provides a 
uniquely accurate measure of the charge carriers injected 
into the system upon K- evaporation, which can be ex- 
pressed in terms of the usual hole doping parameter p, 
i.e. the number of holes per planar copper away from 
half filling (p = 0). 



Electrostatic effects from oxygen and potassium 
doping 

In the following we address the question of the sen- 
sitivity of bilayer splitting to doping and crystal struc- 
ture for various oxygen concentrations in bulk YBCO, as 
well as to the electrostatic effects associated with differ- 
ent 0^~ fillings in the CuO chain layers sandwiching a 
Cu02 bilayer, or the presence of the surface either as- 
cleaved or potassium evaporated. All these configura- 
tions have been investigated by density functional the- 
ory (DFT) calculations using the linearized augmented- 
plane-wave method implemented in the WIEN2K code 
[34 ; exchange and correlation effects were treated within 
the local-density approximation (LDA) [35] . As a further 
check, we have also done calculations using the DFT code 
SIESTA [36 , with the norm-conserving Troullier-Martins 
pseudopotentials [37 and the double-^ singly-polarized 
basis set; the exchange and correlation effects were again 
treated within LDA, after Ceperley and Alder [38 . The 
results obtained with the two methods are in very good 
agreement, with a <0.02eV variation in bilayer splitting; 
the WIEN2K values are summarized in Tab.[ll 

The YBCOe+x band structure calculations, for x = 1 
and 0, were performed using the crystal structures re- 
ported by Jorgensen et al. [39 for x = 0.93 and 0.09, 



respectively. To study the effect of a non-symmetric po- 
tential on bilayer splitting in a direction parallel to the 
c-axis, we calculated the electronic structure of YBC06.5 
using an artificial layered structure consisting of alter- 
nating full and empty chain layers. As one can see in 
Tab.[l[ the bilayer splitting is only weakly dependent on 
doping. The influence on bilayer splitting of the surface 
termination, either of BaO or CuO chain type, can be es- 
timated with YBCO7.0 slab calculations. We adopted an 
ideal bulk crystal structure and periodic slab geometry 
with 18 A separation between the slabs. The slabs are 3 
unit cells thick and stoichiometric, with BaO termination 
on one side and CuO chain termination on the another, 
which allows us to account for the effects associated with 
the polar catastrophe [40]. In the K-deposited YBCO7.0 
slab calculations the K-position is relaxed, whereas the 
YBCO atomic positions are kept fixed according to the 
bulk structure; the computed distance between K and 
BaO (CuO chain) surface is 2.57 A (2.5 A). As shown in 
Tab.jTj also for 'clean' and K-deposited BaO and CuO 
surface terminations we observe a negligible effect on the 
value of bilayer splitting. We can thus conclude that, 
within DFT, the Cu02 bilayer band splitting is almost 
independent of the doping as well as the change in po- 
tential introduced by different filling of the CuO chain 
and the presence of the surface, whether as-cleaved or 
K evaporated. For the purpose of a more direct com- 
parison with the existing literature, in the paper we will 
refer to the values of nodal and antinodal bilayer splitting 
t^^^{k) as obtained for YBCO7.0 in bulk calculations. 



Nodal quasiparticle weight estimate 

Formally, the degree of quasiparticle integrity is re- 
vealed by Zk = jAcoh{k^cjo) duj^ i.e. the integrated spec- 
tral weight of the coherent part of the single-particle spec- 
tral function probed by ARPES [161 Hi!- Experimen- 
tally, while in the optimally-to-overdoped regime Zk is 
believed to be finite - yet quantitatively undetermined - 
at all momenta both above and below Tc, the situation is 
much more controversial in the optimally-to-underdoped 
regime [16l[T7]. In the following, as in the rest of the pa- 
per, we concentrate on the behaviour of the nodal quasi- 
particles because they are better defined over a larger 
portion of the phase diagram. In addition, since they 
are not affected by the opening of pseudo and supercon- 
ducting gaps, they are not characterized by as dramatic 
a temperature dependence as the antinodal quasiparti- 
cles; they can therefore be studied at low temperatures, 
while still providing insights on the normal liquid under- 
lying the emergence of HTSC. To illustrate the limita- 
tions even in the nodal region in the quantitative estimate 
of the quasiparticle integrity from the spectral weight of 
Acoh{k^ oj) as measured by ARPES, we present a detailed 
analysis in Fig.[5| As shown in Fig.|5^, using YBC06.34 
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FIG. 5: Nodal quasiparticle integrity from spectral weight ratios, (a) ARPES spectra at the nodal Fermi wavevector 
kp.N from YBC06.34. (b) Nodal spectral- weight ratio (SWR), defined as /^^ I{kF,N,(jo) duj/ f~^y I{kF,N,(jo) duj, where 
I{k, uj) = Io{k)f{(jj)A{k, (jj) is the intensity measured by ARPES; the data are presented versus doping for progressively increasing 
integration windows (as indicated in c). In (c), the same SWRs are renormalized to the p = 0.37 value, (d) Nodal quasiparticle 
renormalization Zn as determined from the observed bilayer splitting and the SWR estimated at Ep (the latter was rescaled 
so that the two curves match in the 0.23 — 0.37 doping range), together with the mean-field relation 2p/{p-\-l) [131 [14]; thick 
gray and orange lines are guides to the eye. For the bilayer splitting data, the error bars are defined from the B-AB MDC 
fits when bilayer splitting is detected, and from the experimental resolutions otherwise; for the SWR data, they are calculated 
from the spread in the SWR values for integration windows smaller than [0.1, —oo] in binding energy (BE). 



as an example, a clear peak is detected at large dop- 
ing values, which loses intensity upon under doping and 
eventually disappears. The first practical challenge in 
estimating Zj. = jAcoh{k,uj) du is that ARPES does not 
measure directly A{k^uj) = Acoh{k,cjo)^Aincoh{k,uj) but 
rather I{k,uj) = Io{k)f{uj)A{k,uj). This is proportional 
to A{k, uj) via the Fermi-Dirac distribution function f{uj) 
and the one-electron dipole matrix-element Io{k)^ which 
however remains experimentally undetermined (this is 
the so-called 'matrix-element effect' that makes the abso- 
lute value of the photoemission signal a strong and gen- 
erally unknown function of polarization, photon energy, 
and experimental geometry [l6l[T7])- To circumvent this 
problem, since the dipole matrix element is the same for 
both Acoh and Aincohi one can try to evaluate experi- 
mentally Zk = J Icoh{k, cjo) duj/ f I{k, oo) duo. This is pre- 
cisely the strategy followed for instance in Ref.[2T]and[22l, 
which however leads to two additional challenges: (i) As 
already obvious from the data in Fig.|5^, Icoh is strongly 
dependent on doping and progressively more difficult to 
identify upon underdoping; a variety of phenomenologi- 
cal methods have been applied to estimate and remove 
the incoherent background from /(/c,co'), either based on 
the fit of the background by ad-hoc polynomial, linear, 
and integral functions [21] [22] , or on the estimate of the 
background from the spectral weight detected beyond kp 
[24| [4T ] (i.e., from in principle unoccupied regions of mo- 
mentum space), (ii) Aincohi and thus the integral of 
the complete A{k^uj)^ extend over a large energy scale 
given by U /2 c^ 4 eV lU [TTJ [15] ; they cannot be properly 



quantified due to the overlap with other electronic states, 
such as oxygen bands starting already at 1-2 eV BE in all 
cuprates and, in YBCO in particular, also BaO and CuO 
chain bands. As a consequence, even if one could esti- 
mate J Icoh{k,uj) duj, an absolute quantitative value for 
Zk would still be lacking due to j I{k^uo) duo. 

To follow the evolution of the low-energy spectral 
weight versus doping, and attempt an heuristic distinc- 
tion between its coherent and incoherent components, in 
Fig.jsj) we plot the nodal spectral weight ratio (SWR) 

defined as /^^ . I{kF,N^^) duo/ J^ ^y I{kF,N^^) duj, for 
various integration windows extending from above Ep 
(— oo) to a finite binding energy below Ep (BEmin)- Al- 
though the SWR is numerically dependent on the ar- 
bitrary choice of the full integration window ([0.8, —oo] 
in this particular case to avoid the large contribution 
from the oxygen-related features), it is obvious that the 
narrower integration windows will capture more directly 
the degree of quasiparticle coherence. As one can see in 
Fig.[5]3, and even more clearly in Fig.jSJ^ where the SWR 
data have been renormalized to the p = 0.37 value, there 
is a large systematic dependence on the choice of the in- 
tegration window. Although it is not possible to reach 
any quantitative conclusion based on such kind of data, 
a clear qualitative change occurs when the window ex- 
tends beyond the 50 meV BE range: for BEmin ^ 50 meV 
the SWR consistently vanishes below pc^0.15, while it 
becomes finite at all doping - and progressively larger 
- for wider integration windows (eventually approach- 
ing unity once BEmin gets closer to 0.8 eV). Since, as 



shown in Fig.|5^, the sharpest component of the nodal 
spectra is located well within the first 50 meV across the 
whole phase diagram, it is natural to associate the low- 
energy SWR with a relative measure of the quasipar- 
ticle integrity Zn] remarkably, the doping evolution of 
this low-energy SWR is also reminiscent of the quasipar- 
ticle renormalization inferred from the observed bilayer 
splitting through the relation Z 



N 



B,AB /^^LDA 



When the two results are plotted together (after the nec- 
essary rescaling of the SWR, so that the two curves match 
in the 0.23—0.37 doping range), the similarity is striking. 
This validates the quantitative determination of Zn from 
the evolution of bilayer band splitting and confirms our 
conclusion of a rapid loss of nodal quasiparticle integrity 
in the p:^ 0.10 — 0.15 doping range, possibly consistent 
with the results of cluster dynamical mean-field theory 
which predicts a rapid suppression of Zn in this same 
doping range p!5l . 
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